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ABSTRACT 

We present Spitzer IRS and IRAC observations of the young supernova remnant 
E0102 (SNR 1E0102-7219) in the Small Magellanic Cloud. The infrared spectra show- 
strong lines of Ne and O, with the [Ne II] line at 12.8 //m having a large velocity 
dispersion of 2,000-4,500 km s~^ indicative of fast-moving ejecta. Unlike the young 
Galactic SNR Gas A, E0102 lacks emission from Ar and Fe. Diagnostics of the observed 
[Ne III] line pairs imply that [Ne III] emitting ejecta have a low temperature of 650 
K, while [Ne V] line pairs imply that the infrared [Ne V] emitting ejecta have a high 
density of ~10^ cm~^. We have calculated radiative shock models for various velocity 
ranges including the effects of photoionization. The shock model indicates that the 
[Ne V] lines come mainly from the cooling zone, which is hot and dense, whereas [Ne II] 
and [Ne III] come mainly from the photoinization zone, which has a low temperature 
of 400 - 1000 K. We estimate an infrared emitting Ne ejecta mass of 0.04 Mq from the 
infrared observations, and discuss implications for the progenitor mass. The spectra 
also have a dust continuum feature peaking at 18 ^m that coincides spatially with the 
ejecta, providing evidence that dust formed in the expanding ejecta. The 18 ^m-peak 
dust feature is fitted by a mixture of MgSiOa and Si dust grains, while the rest of the 
continuum requires either carbon or AI2O3 grains. We measure the total dust mass 
formed within the ejecta of E0102 to be ~0.014 Mq. The dust mass in E0102 is thus a 
factor of a few smaller than that in Gas A. The composition of the dust is also different, 
showing relatively less silicate and likely no Fe-bearing dust, as is suggested by the 
absence of Fe-emitting ejecta. 

Subject headings: sup ernovae: general - infrared: general - dust:ISM - supernova rem- 
nants:SNR 1E0102-72.3 
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Introduction 



Young supernova remnants (YSNRs) are the most viable astropiiysical laboratories for the 
study of dust formation, nucleosynthesis of heavy elements, cosmic-ray acceleration, and shock 
physics. Of these, dust formation is the least well-studied despite its invocation in efforts to 
explain some of the most important problems in astrophysics, including the isotopi c anomalies of 



heavy elements in meteorites (Clayton 1982), dust formation in the early Universe (iNozawa et al 
20031 ). the measured abundan ces of Galactic co smic rays (Ellison, Drury & Meyer 1997), and the 
interstellar dust budget crisis (jJones et al.lll994l ). 



Meteoritic and astronomical studies show that presolar, cosmic grains condense in the dense, 
warm stellar winds of evolved stars and in the ejecta of supernovae. Mantles of the pre-existing 
dust in molecular clouds are vaporized as the forming stars and planetary systems heat them. A 
small fraction of the dust survives solar system formation without alteration, protected inside as- 
teroids. The most abundant presolar g rains are SiC, nanodiam onds, amorphous silicates, fosterite 
and enstatite, and corundum (AI2O3) (IMessenger et al.ll2006l ). Some isotopic anomalies of heavy 
elements in meteorites have been attributed to the dust that had condensed deep within expand- 
i ng supernovae and some have the r-process composition commonly associated with Type II SNe 



(jClavton et al.lll997l . 12004 . references therein) 



Recent deep sub-mm observations have also shown there to be galaxies and QSOs with very 
large dust masses (> 10® Mq) at z > 5 (Wang et al. 2008; Beelen et al. 2006; Bertoldi et al. 
2003). The timescales for low mass (AGB) stars to release their dust are too long to explain these 
high redshift systems (Morgan & Edmunds 2003). In contrast, supernovae evolved from massive 
stars produce copious amounts of heavy elements and release them on short timescales. Theoretical 
modeling of the conditions in the supernova ejecta indicates that Type-II SNe are sources of dust 



explosion ( 


Deneault. Clavton. k Heeer IboOS: 


Deneault. Clavton. & MeverlboOS: 


Todini & Ferrara 


2001; 


Nozawa et al. 


2003). Yet until verv recentlv. there existed little observational evidence that 



this actually occurs. 

There is now sufficient evidence for dust formation in core-collapse supernovae. For SN 1987A, 
this includes dust emission, dust absor ption and a drop i n line intensities for the refractory elements 



that signals that dust is being formed (iLucv et al 



ejecta of SN 1987A was explored theoretically (jKozasa et al 



1989 



199 l|)- Formati on of dust in the expanding 
19891 . Il99lh . The Type II-P SN1999em 



also showed clear signs of dust formation (Elmhamdi et al. 2003). Detections of CO fundamentals 
with Spitzer obse rvations of SN 2005af and Cas A further provide an indirect evidence of dust 



formation in SNe teotak et al. 



2006 



Eho et al.ll2009l l. 



For the young Galactic SNR, Cas A, ISO observations gave evidence for the association of 
the dust with the ejecta by requiring a mixture of dust grains that are not typical of the ISM, 
while showing strong spatial correlations between the high-velocity infrared emission lines, the dust 
continuum emission, and the optically emitting ejecta (Lagage et al. 1996 and Arendt et al. 1999). 
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Submillimeter observations of Cas A and Kepler with SCUBA (jPunne et alj 120031 : iMorgan et al 
20031 ) suggest the presence of large amounts of cold dust (~ 0. 3 — 2Mf7^ a t 15-20 K), but with 



some controversy related to the presence of foreground material. iKrause et al.l (|2004l ) showed that 
much of the 160 /im emission observed with Multiband Imagin g Photometer for Spitzer (MIPS) is 
foreground material, suggesting there is no cold dust in Cas A. IWilson &: Batrla I (|2005l ). however, 
used CO emission towards the remnant to show that up to about a solar mass of dust could still 
be associated with the ejecta rather than with foreground material. New submm polarimetry 
data confirms that a significant fraction (>30%) of the submm flux originates from within the 
remnant (Dunne et al. 2009). The case for dust formation in supernova ejecta was strengthened 
by Spitzer observations of Cas A, for which line-free dust maps were found to remarkably resemble 
the infrared ejecta- li ne maps of [Ar II], [O VI], and [Ne II]; the estimated dust mass is between 
0.02 and 0.054 Mq (jRho et al.l 120081 ). A similar dust mass is reproduc ed by theoreti cal models 
of dust for mation that include dust de s truction and s tochastic heating (jKozasa I l2009l ) and dust 
destruction (jBianchi &: Schneider 1120071 : iBianchi II2009I ) for Cas A. 



Another excellent source to study dust formation in SNe is the young supernova remnant 
(SNR) IE 0102-7219 (SNR B0102-72.3; E0102 hereafter) in the Small Magellanic Cloud. As 
summarized in Table E it is located at R.A. 0l'"04'^2.57'' and Dec. -72°01'52.3" (J2000), and 



has a diameter of 44", which corres 



fflilditch. Howarth fc Harries 



2005 



jonds to ^ 13.04 pc for a distance of 61.1 kpc to the SMC 



Keller &: Wood 1120061 ). Various wavelength images of E0102 



are shown in Figure [TJ Highly enriched abundances of oxygen and neon relative to carbon and 
magnesium (Blair et al. 2000), as is characteristic of massive progenitors, place E0102 in the class 
of "oxygen-rich" SNRs, along with Cas A. The optical ejecta in E0102 have substan tial velocities 
that exceed 1000 km and imply a dynamic age of about 1000 yr ( Tuohyl Il983l ). The X-ray 
proper motion of the blast wave measured by Hughes et al. (2000) is also consistent with this 
age. More recent HST measurements of the ejecta filament proper motions imply a kinematic age 
of 2050ib600 yr that is a factor of two higher (Finkelstein et al. 2006), but the proper motion 
measurements are consistent within their quoted errors. The compositi on of the optically emitting 
eje cta suggests a inassiv e progenitor that underwent either a Type lb (jBlair et al.ll200d ) or Type 
Ic ([Flanagan et al.ll2004l ) supernova explosion. The optical expansion rates place E0102 currently 
in transition from free expansion to the Sedov phase. The low extinction toward E0102 compared 
with that of Galactic SNRs makes it possible to observe the infrared (IR) emission wit hout contam- 



i natio n by emission from clouds along the line of sight (as seen, for example, in Cas A. lKrause et al 



(|2004l )). The low metallicity environment of the SMC also makes E0102 a useful Type II prototype 
to compare with SNRs at higher redshifts. 

In the infrared, the 24 fim Spitzer MIPS image of E0102 presented by Stanimirovic et al. 
(2006) shows a filled morphology with two prominent elongated filaments that resemble structures 
seen in the X-ray image; the remnant was not detected at 8 or 70 fj,m. The infrared emission 
appears to be mainly associated with reverse shocks of the hot X-ray gas. An IR emission peak 
around 24 /im would suggest the presence of hot dust in E0102 with T~120 K, with most of the 
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dust being centrally located, based on the 24 /im image. Assuming the dust originates in the ISM 
and using typical ISM dust absorption coefficients ([Praine &: Lee II 19841 ). the estimated dust mass is 
8 X 10~^Mq. Even if all the dust were actually formed in the explosion, this dust mass is lower than 
might be expected based on some theoretical models. Most importantly, since the 24 /im emission 
contains both line and continuum, not photometry b ut spectros c opy i s required to estimate the 
mass and composition of dust as was demonstrated by lRho et al.l (|2008l ). 



In this paper, we present Spitzer spectral and imaging observations of E0102 with IRS and 
IRAC. We find that the infrared ejecta emission is dominated by Ne and O and detect a dust 
feature peaking at 18 /um which is spatially coincident with the ejecta emission. A high portion of 
the dust clearly comes from ejecta, rather than being all associated with circumstellar/interstellar 
material. We show that some of dust is at lower temperature and a total dust mass is higher than 
inferred from earlier observations. We discuss physical conditions of infrared emitting Ne and its 
inferred ejecta mass using emission line diagnostics and shock models. 



2. Observations 

We performed an IRS staring observation toward the southeastern shell of E0102 (R.A. 0l'^04™04.04'^ 
and Dec. -72°02'00.5", J2000, see Figure |2]) as a part of our Young SNR Spitzer GO program (PI: 
Rho). The Long Low (LL: 15-40 /im) IRS data were taken on 2005 August 14 with 6 cycles of 30 sec 
exposure time; this yields a total exposure time of 360 sec for the first and second staring positions. 
The Short Low (SL: 5-15 /xm) IRS observations were made with 3 cycles of 60 sec exposure time 
and one cycle covers 2 dither positions; this yields a total exposure time of 360 sec per sky position. 

The IRS spectra were processed using the S15.3 pipeline products as follows. First, we removed 
rogue pixels (bad pixels mainly due to cosmic rays) by subtracting nearby empty sky fields with 
the same exposure time. We took the empty sky data from our observation in the alternate orders 
or from archival data (AORs 14706944, 14707712, 14708480, and 14708824). Any remaining rogue 
pixels were removed by using the SSC tool IRSCLEAkEI. Second, we extracted the spectrum using 
SPICE, and removed fringes. Third, we applied aperture and slit-loss corrections for diffuse emission 
(also see Tappe, Rho, &: Reach 2006). The reduced spectrum is shown in Figure El Additionally 
we subtracted a local background spectrum using the relatively emission-free regions right outside 
the SNR in our slit data (specifically, the region outside the eastern side of E0102 at RA. 16.060, 
Dec. -72.039). The background-subtracted spectrum is shown in Figure HI Comparison of the two 
spectra is described in §3.1. The staring mode yields 4 coadded spectral images (positions 1 and 
2, LLl and LL2) which we merged to make a position- wavelength map as shown in Figure |5l The 
location of E0102 on this map is marked by the appearance of both continuum emission and strong 
[Ne III] (15.55 /xm) and weak [Ne V] (23.3 fim) emission lines. 



http://ssc.spitzer.caltech.edu/archanaly/contributed/ 
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We also used archival IRS mapping data from the IRS Legacy SMC program (PI: Bolatto). 
These maps were made for large regions of the SMC and included E0102. The observations used 
1 cycle of 30 sec exposure time with 5 arcsec overlap, yielding a total of 60 sec per sky position. 
These mapping data were 6 times shallower than our GO IRS staring mode data. 

We made a deep IRAC observation as a part of our GO program. The IRAC observations 
were made for a full array with 30 sec frame time and a 9 position random dither pattern, yielding 
36 sec exposure per sky position. The data were taken on 2005 June 12. IRAC bands cover the 
wavelength ranges 3.2-4.0, 4.0-5.0, 5.0-6.4, and 6.4-9.4 //m, respectively, with a FWHM of 1.7". 



3. Results 

3.1. IRS Spectra: Continuum and Ejecta Lines 

The IRS spectra of E0102 from the southeastern shell in Figure[3]show both line and continuum 
emission. The detected lines include [Ne II] (12.8 //m), [Ne V] (14.3 //m), [Ne III] (15.5 /im), 
[Ne V] (24.3 /im), [O IV] and/or [Fe II] (26 ^m), and [Ne III] (36.07 //m). The local-background 
subtraction effectively removed the [S III] (33.5 /im), [Si II] (34.8 //m), and [S III] (18.7 /xm) lines, 
but this background is likely to be contaminated by the nearby H II region of N76. When we 
generated a [Si II] map from the IRS cube, E0102 appeared as a [Si II] emitter. The lack of 
clean local background data thus makes it unclear whether there are actually any S and Si lines 
associated with E0102. The position-wavelength LL maps in Figure [5] show that the SNR E0102 
emits continuum emission, [Ne V], [Ne III], and [O IV] (or [Fe II]). The 26 /im blend is likely to be 
dominated by [O IV], and unlikely to contain any [Fe II], because other [Fe II] lines such as those 
at 17.9, 6.7, and 5.3 nm are not detected. The infrared spectrum of E0102 is thus dominated by 
Ne and O lines. 

The detected lines and their brightnesses are summarized in Table [2l We measured the line 
brightnesses using two methods: a Gaussian fitting, and integration of line fluxes over wavelength 
ranges. The latter is required in order to account for red- and blue-shifted velocity dispersions. 



To o btain unreddened line brightness, we used an extinction value (Av) of 0.08 mag (|Blair et al 
19891 ). which is equivalent to 2 ( <8)x 10^*^ cm~^. This value is consistent with those of X-ray 



measurements (ISasaki et al.ll2006l ). The extinction corrected values of A^ are also shown in Table 
[TJ Because the extinction is very small, the unreddened line brightnesses differ by less than 0.3%. 

The E0102 spectrum shows a dust feature between 16 and 20 fim peaking at 18.5 ^um (we 
call it the 18 /im-peak dust feature, hereafter). This feature occurs in the wavelength range of a 
prominent silicate feature, but it is very broad and extends to shorter wavelengths. We present a 
multi-composition fit in §4.1. 
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3.2. IRAC images 

Faint emission, particularly from the southeastern shell, is seen in our IRAC 8 fim image of 
E0102; the 8 /xm emission corresponds well with the MIPS 24 //m emission, as shown in Figure 
[6l The flux range in the image is 0.15-0.2 MJy sr~^. However, there is also some 8 /um polycyclic 
aromatic hydrocarbons (PAH) emission from the nearby H II region of N76 to the southwestern 
side of the SNR and a large patch of diffuse emission outside the southwestern shell of the SNR. We 
did not detect 8 iim emission with IRS, but IRS has poorer sensitivity than IRAC, with a limit of 
0.24 MJy sr-i for our IRS staring observation compared to 0.15 MJy sr ^ for IRAC. Because the 
IRAC flux is measured relative to the nearby background, we believe that the IRAC 8 fim emission 
detected towards E0102 comes from the SNR. Moreover, dust emission is expected in the ejecta 
and any emission from E0102 at 8 fim is likely to be continuum emission. A deep IRAC image, or 
a more sensitive observation with future missions may be helpful in unambiguously determining if 
E0102 does emit at 8 fiui. No emission of E0102 at 3.6, 4.5 and 5.8 ^m is detected. 



3.3. Line and Dust Maps 

We generated line and continuum maps using the IRS mapping data, which covered the SNR 
E0102 and the nearby H II region N76. The maps are shown in Figure [71 The line maps were gen- 
erated for [Ne II] , [Ne III] , [O IV] and [S III] by subtracting the continuum baselines at neighboring 
wavelengths. Infrared ejecta emission from [Ne II], [Ne III], [O IV] and [Si II] (Figure [7]) coincides 
with the optical [O HI] ejecta emission as well as with X-ray ejecta emission. In spite of the lower 
angular resolution for the infrared images, this can be seen clearly in the overlay of the [ Ne II] image 



with the optical [O HI] contours (from Blair et al. 1989) and the X-ray Chandra image (jCaetz et al 



2000) in Figure [8] . The infrared emission is strong at the southeastern shell and center, both of 
which coincide with ejecta emission seen in X-ray and optical images. The infrared line emission is 
weaker in the northern and northwestern regions, which are, however, also associated with X-ray 
emitting ejecta. The ejecta line profile of the northwestern region shows significant line-broadening 
(see §3.4 for details). The infrared continuum image is most similar to the infrared [O lV]-|-[Fe II] 
(at 26 fj,m) line image, except in the northwestern shell. However, this continuum emission at the 
northwestern shell nicely follows the X-ray ejecta shell as shown in Figure [TJi, and is located inside 
the forward shocked material traced by the radio shell (see Figured]). The dust emission is thus 
seen to be spatially correlated with the ejecta emission. 

We compared the line and continuum maps of E0102 with those of the H II region, N76, to 
its southeast (see Figure [9]) , in order to examine their infrared colors and to assess the effect of 
background contamination by N76. E0102 shows stronger [Ne II] and weaker [O IV] emission than 
N76. The Ne lines of the SNR are thus confidently measured. The continuum and [S III] emission 
of N76 is strongest at its outer edge, while [O IV] is strongest at the inner shell of N76, as shown 
in Figure m Since E0102 has its strongest [Ne II], [O IV], and dust continuum emission in the 
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southern and eastern parts, we are confident that the Ne, O and continuum measurements are not 
contaminated by N76. The [S III] hue image shows that the SNR also emits [S III] emission as well 
as the H II region. Therefore, the [S III] line in the IRS spectra likely belongs to the SNR. We use 
the spectrum without local-background subtraction in Figured] for the spectrum of E0102, but the 
difference from the local-background subtracted spectrum is used in estimating the uncertainty of 
the spectrum. 

Infrared ejecta seem to most closely correlate with optical emission. Figure [8] demonstrates the 
correlation among the X-ray, optical and three of our infrared ejecta maps. Infrared [O IV] is a fair 
direct comparison to the optical [O III] line maps, but the [O IV] image has a factor of few worse 
angular resolution than that of the [Ne II] image. [Ne II] and [O IV] correlate with optical [O III] 
images in the southwestern shell and center. It is likely the infrared ejecta are denser, cooler ejecta 
than the optical ejecta. In the future, the infrared and optical data should be combined to more 
accurately estimate of the ejecta mass, because a wide range of ionization states are observed in the 
remnant. However, it will be more meaningful to combine them within similar resolution images, 
since the eject a knots have very s mall scales not only in optical but also in infrared as shown in the 



case of Cas A (jEnnis et al.ll2006l ). The dust continuum emission correlates with [O III] and X-ray 
emission for most of E0102 as well as with infrared ejecta in the southwestern and central parts of 
the SNR. This may indicate that the freshly formed dust is correlated not only with infrared ejecta, 
but also with optical and X-ray ejecta (recall that the infrared continuum shell is still located inside 
the forward shock so it is less likely due to CSM/ISM material). Dust forms in the ejecta when 
the gas cools below a temperature of 1000-2000 K. After the reverse shock encountered the ejecta, 
it sufficiently heated both gas and dust, which emit the infrared emission lines and continuum, 
respectively. 



3.4. Doppler-shifted lines and maps 

The measured [Ne II] line width of the staring-mode spectrum (SE shell position in Figure [2]) is 
0.1529ib0.0023 //m (Table [3]). To examine the velocity broadening, we extracted SLl spectra from 
three additional positions in E0102-center, east, and Scenter (see the boxes of Figure [T0|) . The 
three spectra show different line widths, as shown in Figure [TTl The doppler-shifted line maps in 
Figure [TOl with the locations of the three spectra overlaid show that the central position has both 
blue- and red-shifted emission, while the eastern position has blue-shifted emission. The properties 
of the line widths at these positions are summarized in Table [3l At the center, the measured ~0.12 
/im line width is comparable to those measured for the [Ne II] calibration sources G333.9, NGC6720 
and NGC7293. By contrast, the widths in the center and E spectra do indicate line-broadening. 

Fitting one Gaussian component, the line widths correspond to velocity dispersions of 5226ib304, 
3579±54 km s~^, and 3433ib64 km s~^ for the center, SE shell, and E, respectively. After taking 
into account the spectral resolution, we can estimate the true velocity dispersion for each position 
(Table [3]). The center, E and SE shell positions show the true velocity dispersion of 4446]*^ gg'^. 
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2060ljQg, and 2295lg| km s ^, respectively. The observed velocity shifts and dispersions are com- 
parab le to those from optical measurements from -2600 to 3640 km s~^ (jTuohvl 119831 : iBlair et al 



2000) • The broad line profile for the center can also be fitted with two Gaussian components. Using 



a rest wavelength of [Ne II] of 12.8135 fim, the ejecta material must then be red-shifted by 2615 
km s~^ and blue-shifted by 641 km s^^. A confirmation with higher-spectral resolution data is 
desirable because of the limited spectral resolution of IRS low-resolution spectra. 



3.5. Ne Line-Flux Ratio 



In the line list in Table [2l the ions [Ne III] and [Ne V] offer pairs of lines which are suitable 
for line diagnostics in order to constrain densities and temperatures. The measured line flux ratios 
for [Ne V] A 14.3/24.3/im are 1.76±0.11, and for [Ne III] A 15.6/36 ^m 54±28. The ratio 15.6/36 
^m is assigned an additional 50% systematic error for the 36 fj,m line intensity, because this line 
falls on the degrading part of the array. To constrain temperatures and densities, we calculate the 
line intensities and ratios of [Ne V] and [Ne III]. We solve the excitation-rate equations including 
collisional and radiative processes as a matri x using 5 energy le v els fo r [Ne V] and 3 levels for [Ne 
III]. The input atomic data were taken from Griffin &: Badnell (j200ol and include temperature- 
dependent collisional strengths. The solutions of the [Ne V] line diagnostics d iffer between the 
newe r temperature-ind ependent collisional strengths of iGriffin &: Badnelll (|2000l ) and older values 
(e.g.. lOsterbrock Ill989l ). while those of [Ne III] are almost the same. The line diagnostics are shown 
in Figure [T2j The density and temperature jointly obtained from the ratios of [Ne V] and [Ne III] 
are 7x10^ (>9000) cm~^ and 612^200 K, when assuming that both [Ne III] and [Ne V] gas come 
from the sar ne gas. Such a high density seems to be characteristic of SN ejecta as observed in 
optical lines: IChevalier &: Kirshner I (|l979l . Il978l ). for example, suggested that the optical [S II] line 



in Gas A originated from ejecta with a density of 10^ cm 



4. Discussion 

4.1. Dust Spectral Fitting and Dust Mass of E0102 

The IRS spectra of E0102 have a prominent dust feature peaking at 18 fim which coincides 
spatially with the emission from the infrared-emitting ejecta. This is clear evidence that dust is 
forming in the ejecta. To determine the dust composition and mass, we performed spectral fitting 
to the IRS dust continuum using the deep IRS spectra shown in Figure [3l Note that the wavelength 
range used for the fits is 5-40 fim, and that only the continuum emission is used to estimate the 
dust mass. 



http://www-cfadc.phy.ornl.gov/data_and_codes/aurost/aurost_excit/ 



- 9 - 



First, we estimated the contribution of the synchrotron emission to the infrared spectrum 
using the radio fluxes an d spectral index. The radio flux at 408 MHz is 0.65 Jy for the entire 
SNR (|Amv &: Ball Ill993l ). We measur ed the radio surface b rightness at the IRS position to be 0.58 
MJy sr~^ from the 408 MHz image of lAmv &: Ball I (|l993l ). For a radio spectral index a = —0.70 
(jAmv Sz Ball Ill993l ). where S,^ oc i^", the expected synchrotron fluxes are 1.4 x 10"'*, 3.7 x 10~^, 
and 5.5 x 10"^ MJy sr'^ at 5, 20 and 35 /zm, respectively. This is a very small contribution to the 
infrared continuum (<1%). 

The dust continuum was then modeled with the Planck function i?j/(T) multiplied by the 
absorption efficiency {Qabs) for various dust compositions, with the amplitude and temperature 
of different components allowed to vary. The fitting technique and mass estimation method is 
the same as that was used for Cas A (Rho et al. 2008). The spectrum shortward of 10 //m is 
highly uncertain due to background subtraction and instrumental noise. Three dust components 
are needed to fit the IRS spectrum of E0102. 

The 18 /im-peak dust feature includes silicates, but it cannot be reproduced by the combination 
of C, Mg-rich silicates, and AI2O3 grains that are most likely to be associated with the Ne- and 
0-rich ejecta. A component with a feature around 16 /xm is also needed. Candidates includ e MgO 
and Si. Whereas the MgO feature calculated using the optical constants of Hofmeister et al. (j2003l ) 
turns o ut to be too narrow, amorphous Si, which can condense in the C-rich layer in the primordial 
SNe II (INozawa et al.ll2003l ). can successfully reproduce the observed feature. 



The rest of the continuum requires either carbon (Model A) or AI2O3 (Model B) dust, as 
shown in Figures O and [T5j We favor carbon or AI2O3 over solid Fe dust, because we expect 
carbon or AI2O3 dust to be present where the Ne and O eiecta lines are dominant: Ne, Mg, and Al 
are all carbon-burning nucleosynthesis products (jWooslevlll995l ). MgSi03 and AI2O3 are condensed 
in the ejecta of Ne, Mg, and Al, which are also carbon-burning products; the O and Al are found 
in the outer layers of ejecta, which is where Si and carbon dust grains condense. The dust and 
electa compositi ons seen in E0102 are similar to those of the weak 21 /im dust layer in Cas A 
(jRho et al.ll2008l ). E0102, however, does not have the same silicate dust features seen in Cas A. 
Also, unlike Cas A, the dust composition in different ejecta layers in E0102 cannot be spatially 
resolved, at least partly because observations of E0102 suffer from limited angular resolution due 
to its 18-times greater distance. 

The dust compositions, temperatures, and masses obtained from the best fits are summarized in 
Table [H The estimated dust masses within the IRS observed slit are 5.0x10^3 and 5.5 xlO"^ M© for 
Models A (carbon) and B (AI2O3), respectively. To determine the sensitivity to dust temperature, 
in Model C we performed the fit allowing the temperature of AI2O3 to vary independently. Model 
C produced a low temperature (< 45 K) that cannot be constrained by our spectra, which lack 
long-wavelength data; therefore, we don't consider Model C further and fixed the temperature of 
AI2O3 in Model B to that of carbon in Model A. To explore the effect of including many minerals, 
in Model D we performed a fit including all the most common dust compositions (i.e., MgSiOs, 
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Si, AI2O3, carbon, and Mg2Si04; see e.g.TOl, N03). The mass is relatively well constrained by the 
observations (within 20-40%). The fit of Model D is shown in Figure [T6l However, the reduced 
~ 2 is still relatively high, largely due to the residuals near where the gas lines were subtracted. 
(Note that some of lines show kinematic broadening.) 

To estimate the total dust mass, we correct the observed mass (within the high-sensitivity 
spectral slit) for the fraction of the SNR covered by the IRS slit using a continuum map (between 
21 and 23.5 fj,m) generated from the lower-sensitivity spectral cube. This approach is justified 
because the spectra from the cube did not show significant variation across the SNR. The correction 
corresponding to the total flux of the entire SNR relative to the flux of the regions covered by the 
IRS slit is a factor of 2.9. The total dust masses for the entire SNR are then 0.015, 0.007, and 0.014 
M.Q for Models A, B, and D, respectively. We favor Models A or D over Models B or C, because 
the temperature of AI2O3 could not be constrained. 

MgSiOa, A I2O3 and carbon a re ma j or dust specie s pred icted to be produced in the ejecta 
of SNe by both iTodini fc Ferraral (|200ll ): iNozawa et al.l (12003). Note that presolar AI2O3 grains 
from meteorites have be en inferred to be ani ong the most abundant isotopically-enriched materials 
ejected by Type II S Ne (IClavtqn et al.ll2004l '). In term s of amorphous Mg2 Si04 versus MgSiOs, the 



young SNRs Cas A (IRho et al 



200. 



and SN 2006 



Nozawa et al. 



c (INozawa et al.l 120081') also s how a larger mass 



20031) and iTodini fc Ferraral (l200lh predict the 



of MgSiOa than Mg2Si04, whereas 
reverse. This may suggest that the chemical network favors MgSiOs over Mg2Si04. 



4.2. Ionization States of Ne 

Column densities for [Ne III] and [Ne V] can be determined from the excitation matrix given 

by: 

I{t,j) = (^-^^ NiA{i,j)huij (1) 

where is a line intensity between levels i and j, A'^j is a column density in the level i, A{i,j) 

is spontaneous radiation transition rate, Aj is a column density in the level i, and Uij is the energy 
difference between i and j levels. For level i, the matix to solve the column density Nj is given by: 

T.j^i Ni Aij + T.jNin 'jij = T.j^i Nj Aji + T,j Nj n jji (2) 

44 = ^ exp-'^-^/^^ (3) 

where n is density, T is temperature and 7(i, j) and 7(j, i) the collisional transition and deexcitation 
rates. The total column density of N = SjAj is obtained by solving the matrix. For example, the 
total column density of [Ne V] is the sum of all five levels as shown in Figure [T3l For a temperature 
of 612 K and a density of TxlO'' cm~^, the column densities for [Ne III] and [Ne V] are 6.8x10^^ 
and 5.9x10^^ cm~^, respectively. We estimated the column density for a few sets of densities and 
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temperatures implied by the observed [Ne III] ratio as shown in Figure [T2\ and found that the 
column density changes rapidly depending on the electron density as listed in Table [5j We also 
solved the excitation-rate equation using 2 energy levels for Ne II. The observed line brightness of 
[Ne II] can be used to estimate the ion column density to be 2.7x10^^ cm~^. The column densities, 
ionization potentials, and critical densities of Ne are summarized in Table [5j 

The nucleosynthetic yields of Ne trace the progenitor mass. In order to obtain the Ne ejecta 
mass, we must estimate the column densities of the remaining ionization stages of Ne. First we 
obtain the column densities when we assume local thermodynamic equilibrium, which is described 
below in this section. Later we use shock models to estimate the column densities of ionization 
states as described in §4.3. 

The Saha equation is given by: 



In ( )+36.11413+-lnr-11602.81^ 

(4) 

where N is the number density of atoms in the r (or r+l)th stage of ionization, and U is the partition 
function depending on the temperature, is the mass of the proton, T is the temperature in 
Kelvin, x (lev) is the ionization potential for r to (r+1) (in units of eV). 

We estimated the column densities of [Ne I] and [Ne IV] using the Saha equation repeatedly 
as follows. First, we obtained a gas temperature of 5667 K by using the column densities of [Ne II] 
and [Ne III]. While this is a much higher temperature than obtained from the [Ne III] and [Ne V] 
line diagnostics, it is consistent with the error range for the line ratios (see Figure [T2]) . Second, we 
estimated the neutral [Ne I] column density by using the column density of [Ne II] and a temperature 
of 5667 K. Third, we estimated the [Ne V] column density by using the temperature of 5667 K 
and the Saha equation, and found that the inferred [Ne V] column density was much smaller than 
the observed value of 5.8x10^^ cm~-^. This indicates that most of the observed column density of 
[Ne V] comes from a gas with a temperature different from 5667 K. An additional component with 
a temperature of 9100 K is required to fit the observed [Ne III] and [Ne V] column densities. The 
[Ne IV] column density is then obtained by interpolating between [Ne III] and [Ne V]. Figure [17] 
shows the column density for each ion of Ne from the two-temperature components derived here. 
Column densities for the remaining, more ionized Ne can be calculated for these two temperatures, 
and can be seen to fall rapidly with increasing ionization. We also estimated the column density of 
[Ne III] for (ue, T) = (10"^ cm'^ 350 K) and (10^ cm'^, 215 K). We find that the column density 
varies depending on density, as summarized in Table [5] and Figure [HI and that it does not vary 
much with temperature. A lower electron density yields a higher column density. 
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4.3. Shock Model and Neon Mass 



We have calculated radiative shock models for shocks in the veloc i ty ran ge of 50 - 500 km s ^. 
These models use many of the same assumptions used by iBlair et al.l (j2000l ) . In particular we use 
the same elemental abundances (H = 12.00, O = 16.00, Ne = 15.50, C = 14.50, Mg = 13.90) and 
follow their approach of using a single ra m pressure, npf ^ = 16cm~^100kms~^, for each shock 
speed in the standard case. We also f ollow iBlair et al.l (l200d ) in using a magnetic field strength of 
Bq = A^/uq fiG for our standard case. iBlair et al.l (|200d ) use the results of a range of shock models 
with shock velocities ranging from 30 — 500 km s~^, weighting the emission by the fraction of the 
shock area that is assigned (by assumption) to each shock speed. We take the simpler approach of 
finding the best fitting single shock speed. In addition to the standard cases, we have calculated 
models with ram pressures that are larger by a factor of 1.5, 2 and 3 and with smaller magnetic field 
strengths, down by a factor of 2. The model calculations assume 1-D (plane parallel) steady flow 
and include non-equilibrium ionization, radiative cooling and magnetic pressure, B oc n. We also 
include photoionization from radiation generated in upstream gas and the accompanying heating. 
Since the region of strong ionizing radiation emission is mainly in the hot, ionized portion of the 
shock, while the region of substantial photoionization is downstream, this approach is satisfactory. 
In fact photoionization and heating are negligable within the X-ray/ EUV emitting zones where 



T > 10^ K. These calculations make use of a revised version of the [Raymond Sz SmithI ( 19771 ) 



plasma emission code for the calculation of the hot, ionized post-shock zone and the cooling zone 
in which the gas temperature drops sharply to < 1000 K. 

For the calculatio n of the photqioniz ed zone, we use the Cloudy code (version 07.02.02 of the 
code last described by lFerland et al.lll998l ) with the incident radiation field taken to be the ionizing 
radiation field generated in the hot, collisionally ionized and cooling zones (as calculated using the 
R&S code). We show an example of the ionizing radiation field in Figure [T8l Cloudy assumes 
thermal and ionization equilibrium (in contrast to our calculations using the R&S code), but this 
should not be a very bad assumption for the cool (T 400 — 1000 K), dense gas of the photoionized 
zone. The intensity of the ionizing radiation field depends on the geometrical factor, -Rmaxi re lated 
to the ratio of the lateral extent of the shock front to its thickness (see, e.g. lRavmondlll979l ). We 
have generally used iJmax = 1 in our calculations, but found the best agreement for a model with 
-Rmax = 0.7 indicating small shock fronts as might be indicative of shocks propagating into clumps 
elongated in the radial direction. Figure [19] illustrates the Ne ionization and temperature for the 
two parts of the calculation for one model. 

As illustrated in the figure, [Ne V] emission is entirely from the ionization and cooling zone, 
while [Ne II] and [Ne III] emission comes primarily from the low temperature, T 600 K, pho- 
toionization zone. The electron density and temperature derived in Figure [12] from the intersection 
of the [Ne III] and [Ne V] contours implicitly assumes a uniform temperature and density for the 
[Ne III] and [Ne V] emitting gas. Since the emission comes from different regions with different 
temperatures and electron densities, these diagnostics do not apply to emission from the shock. 
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Under the assumption that the [Ne III] and [Ne V] emission comes from the same gas, we 
calculate the column densities of Ne I, Ne II and Ne III to be 7x10^^, 3.5x10^^, 6.5x10^^ cm~^, 
respectively, and the total Ne column density to be 7.9x10^^ cm~^. To use the results of our shock 
models we need to limit the extent of the photoionized zone. We use the total [Ne II] 12.8 //m 
emission to set the depth of the photoionized zone. The total [Ne V] emission is naturally limited 
by the size of the post-shock ionization and cooling zones. We find, however, that we are unable to 
match both the total surface brightness of [Ne V] 14.3 //m emission and the ratio of 14.3 /im/24.3 
/im emission at the same time when the surface brightness is determined as in Table [2j Matching 
the line ratio requires lower shock speeds which in turn leads to more flux. 

This problem can be overcome if the emission is clumped on scales unresolved by Spitzer. Table 
[2] uses the default angular extraction size (10"x20" for LL) to estimate the surface brightness, 
because the SNR structures were not resolved. However, we know from the other young SNRs 
that the angular size of ejecta emitting regions are generally much smaller than the implied ~ 3 
pc. Therefore, w e estimate the angular size (to get the solid angle) using a Ne image of Cas A 



(jRho et al.ll2008l ). From this we estimate the filling factor of the Ne emission over the 20"xl0" 
field of view to be ~ 0.13 in surface area. The [Ne V] 14.3 //m surface brightness of 3.974x10"^ 
erg cm~^ sr~^ thus becomes 3.057(ibl) x 10~^ erg s~^ cm~^ sr~^ after the correction, where the 
errors are estimated using experiments with different surface brightness contours (the filling factor 
ranging from 0.1 to 0.2). 

With the correction for the filling factor we find that we can match both the [Ne V] line ratio 
and fluxes for shock velocities of ~ 200 km s~^, though a better match is achieved for models with a 



higher, ~ a factor of 3, ram pressure than the standard iBlair et al.l (|2000l ) value, as shown in Figure 
[20l Using the filling factor correction and the best fitting model for the [Ne II] and [Ne V] emission, 
we find the column densities of [Ne I], [Ne II] and [Ne III] to be 1.5x10^^, 7.8x10^^, and 1.2x10^^ 
cm~^, respectively, and the total Ne column density of 4.0x10^^ cm~^. The column densities from 
different ionization states of Ne are summarized in Table O 



4.4. Ne Ejecta Mass 

From the column densities for each ion, we can readily calculate the Ne mass by summing over 
all the ions as: 

M = Y^iNi mNo ^ (5) 

where M is the total mass mass, Ni is the column density of ion i, rriNe is the mass of Ne, d 
is the distance and the emitting area. In our calculation, we used a wavelength dependent 0, 
which is equivalent to the line extraction width (the approximate area of the LL extraction region 
is 10"x20"). The mass corresponding to the column densities for each ion is also shown in Figure 
[T71 The total derived Ne mass (within the IRS slit) inferred from excitation and Saha equations 
is 0.22 Mq for ne= 7x10^ cm'^ and T= 612 K, and 0.07 Mq for ne= 1x10^ cm"^ and T= 215 
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K, respectively. For the shock model we find a total mass in Ne of 0.014 Mq which supercedes the 
derived values. 

This is the mass enclosed by the IRS slit only, which is just a portion of the entire SNR. We 
make the correction for the entire SNR using the [Ne III] map at 15.5/um. The total flux is 3.03 
times larger than that covered by the slit, so we obtain a total Ne mass (for all ions including 
neutral Ne) in E0102 of 0.042 M© from the shock model. 

We compared our derived Ne mass of E0102 with that of SN2005af at an epoch 214 days after 
the explosion as a rough consistenc y check. The Ne mass implied from [Ne II] line of SN2005af 
is 2.2x10-3 Mq (|Kotak et ahlbood ). We determine the [Ne II] mass in E0102 to be 3x10"^ M© 
(from the shock model), which is within a factor of 1.5 of that determined for SN2005af, also a 
massive Type II SN. This is a relatively small disparity given that there may be differences in the 
progenitor mass, the respective physical conditions of Ne-emitting gas, and thermodynamic history 
of Ne gas, and also does not account for hotter and more ionized O ejecta that are observed at 
shorter wavelengths. 

We detected no infrared H lines from E0102, though a few do fall in the IRS spectral band. H 
lines are present in th e optical spectra of E0102, but it is not likely that they are associated with 
the supernova ejecta (iBlair et al.ll2000l ). For the purpose of estimating the Ne abundance relative 



to solar (or cosmic) values, however, we adopt the Ha line flux 5.9 4x10 erg s ^ c m ^ and the 



upper limit of II/3 line flux of 1x10^^^ erg s"^ cm"'^ reported bv iBlair et al.l ( 200d ) in order to 
obtain an upper limit on the H column density. At large optical depth, 4 vr j(H/3)/(Np Ng) is a 
consta nt for a given tem perature that varies by only a factor of 5 over the relevant temperature 
range (jOsterbrock Ill989l ). Taking the path length to be the aperture beam size of 1" (0.299 pc), the 
H column N is inferred to be <7xW^^ cm~-^. The Ne abundance in E0102 ([Ne/H]) using the Ne 
column density of 4x10^^ cm^^ is then >1.9xl0~3, while the solar value ([Ne/H]©) is 1.2x10"^, 
so [Ne/H]/ [Ne/H]0 >18 if the optical Ha is associated with E0102. More likely, that is not the 
case, so this is a lower limit for the abundance of Ne relative to solar values. 

We also ca l culate d the column density of [O IV] using 5 levels for [O IV] lines and the atomic 
data of iTaval I (|2006l ) . The column density inferred from the observed [O IV] line brightness 
is 1.5x10^3 cm~^. The equivalent oxygen mass is 8xl0~^ M©. We also estimated the column 
densities and masses of [O I] , [O II] , [O HI] , and [O V] using the derived column density of [O IV] , 
the Saha equation, and assuming that the oxygen gas has the same temperature as the Ne gas. 
The total mass of oxygen in the IRS aperture is thus inferred to be 0.26 Mq. The calculated mass 
is quite sensitive to the second component temperature, however. For example, a temperature of 
7000 K rather than 9100 K reduces the O mass a factor of two. The uncertainty in the estimated 
oxygen mass may thus be a factor of a few. To obtain the total oxygen mass for the entire remnant 
we scale the masses above by the same 3.30 multiplicative aperture correction factor used for Ne. 
The total oxygen mass estimated for the entire remnant is then approximately 0.78 Mq, as shown 
in Table [71 Note that oxygen mass estimate is much less accurate than that of Ne because we only 
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have one observed oxygen line. 



4.5. Ejecta Mass and Nucleosynthesis Yields 



It is evident that the infrared census of the Ne and O ejecta mass is incomplete because the 
ejecta that emit at optical and X-ray wavelengths are more ionized and at higher temperatures. 
Infrared emission traces dense and relatively low temperature ejecta as indicated by our line diag- 
nostics, whereas X-rays trace much hotter and lower density coronal plasmas. The total neon and 
oxygen masses that we estimate from the infr ared observations a nd sh o ck models are compa red to 
the nucleosynthesis yields of massive stars in [Thielemann et al.l (|l996l ). iNomoto et al.l (|l997l ). and 
Chieffi & Limongi (2003) as summarized in Table [71 

The He- and H-like ions that dominate the X-ray emission are infe rred to have masses of 2 
for Ne and 6 Mq for O in the Chandra study of lFlanagan et al.l (|2004l ): these indicate a progenitor 
mass of 32 Mq. These mass estimates, however, ar e subject to uncertainties in the composition 
of the underlying continuum. iFlanagan et al.l ([200J) make their mass estimate a s sumin g that the 
ejecta do not include light elements, but if this is not the case, iFlanagan et al.l (120041 ) note that 
the X-ray derived ejecta mass could be a factor of 20 lower, at 0.1 and 0.3 Mq, for Ne and O, 
respectively. By comparison, our estimates for the Ne and O ejecta masses from t he infrared 



observ ations are 0.042 Mq and 0.78 Mq, respectively. In the optical/UV analysis of iBlair et al 



(|2000l ). relative line strengths and shock models favor a progenitor of 25 Mq, but no attempt is 
made to estimate ejecta masses for the entire remnant. 

If we take the lower set of values for the X-ray emitting ejecta masses, we determine the total 
ejecta mass determined from the combined infrared and X-ray analysis to be >0.142 Mq for Ne 
and >1.08 Mq for O. Taking the higher X-ray values, these masses are 2.04 Mq for Ne and 6.8 
Mq for O. As a comparison, the X-ray emitting Ne mass predicted by our shock model is ~0.84 
Mq, though we note that this model was matched to the infrared Ne line emission only. Taken 
together, the inferred ejecta masses favor a progenitor mass of 30 (25 - 40) Mq as shown in Figure 
[2Tk and Table[71 but are not strongly constraining. The infrared Ne/0 mass (>0.13) ratio by itself 
favors a progenitor mass greater than 25 Mq, as do the total masses of Ne and O when we take the 
lower range of the X-ray emitting ejecta mass corresponding to the presence of light elements in 
the ejecta. This is not surprising given that the infrared ejecta mass makes up a significant fraction 
of the total mass (30% for Ne, and 70% for O) if the X-ray ejecta masses are taken to be at this 
lower limit. Contrariwise, the infrared emitting masses make up a very low fraction of the total 
ejecta mass (2% for Ne and 12% of O) if we take the higher X-ray ejecta masses corresponding to 
pure heavy elements, and the inf erred progenitor mass is then necessarily in agreement with that 
obtained in the X-ray analysis bv IFlanagan et al.l (|2004l ). 



The set of line-emitting ejecta elements is very different for E0102 and Cas A, although both are 
young core-collapse SNRs with optical emission from high-velocity O ejecta. In E0102, the Ar lines 
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are apparently absent in the infrared, and the Si and S hnes are weak. Cas A, by contrast, shows 
very strong Ar emission as well as strong silicate (21 //m-peak) dust emission (see the comparison 
in Table [6l). In addition, Cas A shows unshocked, centrally located Si and S ejecta, whereas E0102 
shows no indication of any unshocked ejecta. The differences in the X-ray line emission from 
these two remnants a r e in line with thos e seen in the infrared. E0102 is dominated by O and Ne 



(jFlanagan et al.l 12004 : ISasaki et al.ll2006l ). while Cas A is dominated by Si, S, Ar, Ca, and Fe. The 
ejecta temperatures in Cas A are higher, and presumably the O is too ionized to emit appreciable 
line emission, although Cas A is also subject to much higher interstellar absorption. And in spite 
of being considerably younger than E0102 (at an age of 330 yr), Cas A has a reverse shock that is 
already traversing the innermost ejecta layers. The key to this difference may be the extensive mass 
loss experienced by the Cas A progenitor, which probably had a binary companion (Young et al. 
2006). Cas A is likely to have exploded at only 4 Mq after having started at ~25 Mq on the main 
sequence (e.g.. Laming & Hwang 2003). This difference in the mass loss history of the progenitors 
would also affect the supernova nucleosynthesis and subsequent evolution of the remnant. 

It is thus possible that the explosion of E0102 simply did not produce a significant mass of 
Ar. Certainly, variations of factors of 5-20 in Ar mass are evident in va rious supernovae (f or 



example, see the comparison of SN 2005af to SN 1987A and SN2004dj in iKotak et al.1 (12006!)). 
Nucleosynthesis models likewise predict a wide range of Ar masses, with more massive progenitors 
producing significantly more Ar. The yields range from 0.001 to 0.07 Mq for 25-40 Mq progenitors, 
with corresponding Ar/Ne mass ratios of 0.001 - 0.1 (Figure [2T|) . If we scale the predicted Ar/Ne 
mass ratios to the total infrared Ne ejecta mass of 6.5x10^'^ Mq, the models predict an infrared- 
emitting Ar mass between 6.5x10"^ to 6.5x10^^ Mq. 

We can use the detection limit of 3x10^*^ erg s^^ cm~^ sr~^ for [Ar II] at 6.985 /um to estimate 
the corresponding mass of Ar in E0102. Using two excitation levels of Ar, the column density of Ar 
is 1x10^" cm~^ if we assume that Ar is at the same temperatures as Ne. Since it is likely that the 
temperature and density could actually be quite different from those values, we consider a range of 
temperatures from 600 to 20,000 K and densities from 10^ cm~^ to 10^ cm~^. Including a factor 
of 3 correction for the IRS aperture, the final upper limit for the Ar mass ranges from 1.5x10^^ to 
1x10"^ Mq. The upper end of this mass range overlaps the range of the expected infrared-emitting 
Ar mass estimated above. It is thus plausible that the current non-detection of Ar is simply be due 
to limited sensitivity. 

We cannot yet exclude the possibility that the infrared-emitting Ar are truly absent, however. 
The physical conditions of the infrared-emitting Ar are highly uncertain. Possibly, the reverse shock 
has simply not reached the deep layers containing Si, S, Ar, and Ca ejecta. It is also possible that 
the Ar might be detected in a more complete set of deep observations. While the IRS mapping 
data are too shallow to detect faint lines, the existing deep IRS staring mode spectrum covered 
only one-third of the remnant. Given that core-collapse remnants often show rather asymmetric 
ejecta distributions, it cannot be ruled out that deep observations at other positions in E0102 might 
reveal Ar ejecta. 
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4.6. [Ne V] lines: Comparison with Other Astronomical Objects 



Among various astronomical objects, highly ionized [Ne V] is not commonly detected in the 
infrared because of its high ionization potential (97.11 eV), whereas lines of both [Ne II] and [Ne III] 
are common and often bright. In galaxies, the detection of [Ne V] 24.3 fim line of comparable 
brightness to [Ne II] is often taken as a signature of active galactic nuclei (jArmus et al.l 120071 ) 
because [Ne V] cannot be readily produced by O stars. In the Galaxy, planetary nebulae emit [Ne V] 



lines which are produc ed by white dwarfs with high temperatures (jRubin 



2004 



van Hoof et al 



200d : lRubin et allllQQTI ). E0102 is the first SNR which is identified as an infrared [Ne V] emitter. 



Here, we compare the physical conditions of ionic line emitting material such as [Ne V] in the 
ejecta of E0102 with those of other astronomica l objects. The [Ne V] 14.3/24.3/im ratio of 1.76 is 
similar to many planetary nebulae (IRubin II2004I ). an example of which is NGC 7 027. E0102 falls at 



the extreme high end of the range (0.2 to 1.8) observed in active galactic nuclei (jPudik et al. 



2007 



Sturm et alll2002l). Both our (see Figure [T2]) and Dudik et al.'s calculations use the same collisional 
strengths from iGriffin &: Badnelll (1200011 . and we are able to reproduce the [Ne V] 14.3/24.3 fim 
ratio plot (Figure 1 of iDudik et al.l 120071 ) . Our measurements of [Ne V] ratios are not affected by 
the extinction correction, because the intrinsic extinction correction is very low at less than 0.1% 
for E0102. In the SNR E0102, [Ne V] lines come largely from collisionally ionized and radiatively 
cooling zones, whereas [Ne II] is dominantly from the photoionized zone. In active galactic nuclei, 
both [Ne V] and [Ne II] lines are shown to be from photoionization. The ratio of [Ne V]/[Ne II] is 
only 0.06 for E0102 and about 1 for AGN. 



4.7. Dust Mass and Implication of Supernova Dust in Early Universe 

Our dust masses for E0102 are 8-18 times greater than the 8 x IO^^Mq mass estimated by 
Stanimirovic et al. (2006) from MIPS imaging, but the latter assumes a temperature of 120 K and 
typical interstellar dust composition. The inferred dust mass will depend on the composition and 
temperature of the dust, which can only be accurately identified with spectral data such as the 
IRS data used here. The difference between our dust mass for E0102 and that of Stanimirovic et 
al. (2006, hereafter S06) thus arises from a number of factors. First, the IRS spectra indicate that 
there are three composition and temperature components. While two of the three components have 
temperatures comparable to the 120 K assumed by S06, the third component had a much lower 
temperature (55-60 K), and as shown in Table HJ dominates the mass. Moreover, because of the 
differences in dust composition assumed, there are differences in the absorption coefficients used 
and also in the inferred density; the latter are generally higher for the grains we used here, with 
MgSiOa, and AI2O3 having solid-state densities of 3.2 and 2.45 g cm~^, respectively, compared to 
more typical ISM grains, e.g., graphite with p = 2.25 g cm^^. 

Our estimated dust mass is also higher than those generally inferred for supernovae. Possibly, 
it is easier to detect dust in young SNRs than in extragalactic SNe simply because the dust in 
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SNe is too cold to be detected in the infrared by Spitzer, whereas the freshly synthesized dust in 
young SNRs has been heated by the reverse shock to observable temperatures. Another limitation 
for estimating the dust mass of extragalactic SNe is that the spatial resolution currently available 
with IR telescopes for observing supernovae is still poor, and that makes it difficult to distinguish 
between fresh dust formed in the ejecta and ISM/CSM dust. In the case of E0102, we were able to 
resolve and separate the ejecta and the forward shocked materials, and found no significant infrared 
emission from the ISM heated dust. 

The estimated total dust mass of 0.007 to 0.015 Mq for E0102 is a factor of 3 lower than the 



mass (from 0.02 to 0.054 Mq) of Cas A (jRho et al.ll2008l ). though of the same order of magnitude. 
We discuss a few possibilities for the discrepancy between our estiinated dust mass from E0 102 



and the dust mass per SN predicted by models (e.g. iTodini Sz Ferrara 



2001 



Nozawa et al.ll2003l ) to 



account for the dust in high red-shifted galaxies. First, there may be more mass in dust at lower 
temperatures. Our wavelength coverage is only up to 40 ^m, and the MIPS 70 ^m emission has not 
been detected due to confusion w ith the nearby H II region of N76 and the limited spatial resolution 
of the 70/im image, as shown by IStanimirovic et al.l (j2006l ). We could not constrain a temperature 
of 41 K with our IRS data, as shown in Model C of TableHl Even if colder dust (<<55 K) exists, our 
observations would not be sensitive to it. The colder dust components have been detected in Cas 
A and Kepler (Dunne et al. 2009; Gomez et al. in preparation) suggests that they could be present 
in other SNRs. Second, there are large uncertainties in dust destruction rates. The dust mass of 
Cas A may explain the lower limit on the dust masses in high redshift galaxies , when assuming 
conservative dust destruction rates (|Rho et al.l l2008l : iMorgan &: Edmundsl 120031 ) . When a more 
moderate destruction rate is t a.ken account, a much higher du s t mass per SN is needed to explai n 
the dust in the early Universe (jDwek. Galliano. &: Jones 1120071 : iGalliano. Dwek. &: Chanial II2008I ). 
The grain disruption which occur s primarily in SN -shocks is known to be too efficient to explain 
the amount of dust in the ISM (jjones et al.lll994l ). However, new infrared observations of SN 
1987A suggest that the d estruction rate of the silicate grains is a factor of 2 smaller than previous 
values (iDwek et al.ll2008l). and that micro-sized grains can survive the passage of a SN shock wave 
(|Slavin. Jones. &: Tielens II2004I ). Third, an inaccurate estimation of progenitor masses from ejecta 
could cause the discrepancy, since the IMF of massive stars is rather steep. As discussed in §4.5, 
progenitor mass depends on accurate estimates of ejecta masses over a broad range of wavelengths. 
Fourth, if E0102 contains Fe-bearing dust, the dust mass would be higher because the density of 
Fe dust {p = 7.95) is three or more times higher than other grains. Since Fe dust (due to the shape 
of Qafes) is featureless like carbon and similar to that of AI2O3 at low temperatures (< 100 K), it is 
difficult to distinguish between Fe, C and AI2O3. While we assume no Fe dust because Fe gas lines 
are not detected, a small contribution of Fe could significantly increase the dust mass. Fifth, some 
grain types may be larger than assumed, which would increase the inferred mass. Sixth, smaller- 
sized dust in t he ej ecta is preferentially destroyed behind the reverse shock as demonstrated by 
Nozawa et al.l (j2007l ): in the case of E0102, this has been ongoing for 1000 yrs. Deep observations 
with future telescopes, in particular, far-infrared/submm observations and theoretical models of 
the time evolution of dust that is formed in the ejecta would be needed to settle this issue. 
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5. Summary 

Using Spitzer IRS and IRAC, we have carried out a study of the infrared emission from dust 
and ejecta in E0102. Our results may be summarized as follows: 

1. We have detected infrared emission from E0102 in both continuum and lines. These include 
bright lines of [Ne II] (12.8 /Lim), [Ne III] (15.5 //m), [O IV] and/or [Fe II] (26 //m) [O IV] and/or 
[Fe II] (26 fim), weak lines of [Ne V] (14.3 fim), [Ne V] (24.3 urn), and [Ne III] (36 /xm), and possibly 
[SII], [SIII] (33.5 urn), [Sill] (34.89 fim). The ejecta in E0102 are dominated by Ne and O in the 
infrared, as in the X-ray and optical. 

2. The [Ne II] lines are broadened by 20 -100 % above the instrumental resolution, implying a 

true velocity dispersion of about 2000-4500 km s~^. This is comparable to velocities measured for 
the optically emitting ejecta, and indicates that the infrared [Ne II] emission is also from ejecta. 

3. The density of the infrared-emitting Ne ejecta is high. The collisionally excited line di- 
agnostics using [Nc V] and [Ne III] line ratios imply that [Nc III] lines come from a low (^^600 
K) temperature and [Nc V] lines from high density region. A shock model with a best velocity 
parameter of 200 km s^^ indicates that the [Ne V] lines mainly come from the cooling zone, which 
is hotter dense region, whereas [Ne II] and [Ne III] lines are mainly from the photoinization zone 
with a cold temperature of 400 - 1000 K. 

4. We detect a broad dust feature at 18 /xm, that we attribute to MgSiOa and Si; the remainder 
of the dust continuum can be fitted with either carbon or AI2O3 grains. The spatial correspondence 
of the dust continuum and the ejecta lines indicates that the dust is freshly formed in the ejecta. 

5. The continuum show strong correlation with ejecta of E0102, showing that dust has formed 
in the ejecta. We estimate a dust mass of 7 x 10~^ to 0.015 Mq that is higher than previous 
estimates by an order of magnitude, but smaller than that of Cas A. 

6. We estimate the mass of infrared-emitting Ne and O ejecta. We estimate infrared emitting 
Ne ejecta mass of 0.014 Mq and a total Ne mass of 1.26 Mq derived from the best shock models. 
Together with the X-ray masses, these imply a progenitor mass of ^^30 M©. 

For the future, detailed modeling to obtain accurate ejecta masses for young remnants such as 
Cas A and E0102 will be helpful to understand the production of elements such as Ar deep in the 
ejecta and their physical conditions in the remnant. More broadly, the goal is to understand the 
progenitors and the explosions that produced these remnants. Infrared observations with Spitzer 
and ISO have allowed us to identify dust features in SNRs and SNe including Cas A and E0102. 
However, the presence of cold dust, which can be only identified in the far-infrared and submm, is 
still uncertain due to limited angular resolution at these wavelengths. In the near-future, we expect 
Herschel observations will significantly advance our knowledge of cold dust. Sensitive observations 
will be very helpful to find and study the ejecta, so deep infrared mid- and near-infrared observations 
with Spitzer and with future missions such as JWST will significantly advance our understanding 
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of the questions of nucleosynthesis and dust formation in supernova explosions. 
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Table 1: Summary of E0102 Properties" 



distance 61.1 kpc (l'~17.7801 pc) 

angular size (diameter) 44" 

physical size 13.04 pc 

shock radius 6.52 pc (22") 

reverse shock radius 4.45 pc (15") 

age 1000 yr, 2050±600 yr 

extinction Ay = 0.08 mag, Ni^=2xl0^° cm~^ 



References are given in the text. 



Table 2. Observed Spectral Line Brightnesses 



Wavelength 


Line 


FWHM 


Surface Brightness" 


Surface Brightness*' 


Aa 






(/im) 


(erg s"-*- cm~^ sr"-"-) 


(erg s"-"- cm~^ sr"-*-) 




12.8145± 0.0010 


[Ne II] 


0.1529± 0.0023 


5.918E-05± 8.774E-07 


5.525E-05 


0.0022 


14.3475± 0.0023 


[Ne V] 


0.1797± 0.0053 


3.974E-06± 1.834E-07 


3.633E-06 


0.0016 


15.5665± 0.0014 


[Ne III] 


0.2027± 0.0032 


2.354E-05± 1.079E-06 


2.491E-05 


0.0019 


24.3410± 0.0059 


[Ne V] 


0.3395± 0.0142 


2.260E-06± 9.676E-08 


2.364E-06 


0.0016 


25.9132± 0.0006 


[O IV] + [Fe II] 


0.3430± 0.0015 


2.077E-05± 9.442E-08 


2.116E-05 


0.0015 


33.5160± 0.0010 


[S III] 


0.3489±0.0040 


2.394E-06±3.023E-08 


2.490E-06 


0.0011 


34.8901± 0.0010 


[Si II] 


0.3459±0.0027 


3.429E-06± 2.920E-08 


3.612E-06 


0.0010 


36.0783± 0.0041 


[Ne III] 


0.2565±0.0197 


4.329E-07±3.595E-08 


5.206E-07 


0.0009 



"Estimated using a Gaussian fitting. 

''Estimated using integration of line fluxes over wavelength ranges. 



Table 3: Properties of Ne Line-broadening 



Location 


RA. 


Dec. 


color 


Wavelength (/im) 


FWHM (pm) 


Velocity (km s ^) 


Velocity^V^, (km s^i) 


Center 


01:04:02.05 


-72:01:53.2 


blue 


12.8081± 0.0057 


0.2235± 0.0134 


5226±304 


44461^53 








(1) 


12.7861±0.0031 


0.16733±0.0070 


3909±164 [641°] 








(2) 


12.9252±0.0045 


0.09391±0.0068 


2196±1640 [2615"] 




E 


01:04:05.01 


-72:01:55.6 


green 


12.7931± 0.0012 


0.1465± 0.0028 


3433±64 




Scenter 


01:04:02.5 


-72:01:59.6 


red 


12.8104± 0.0026 


0.1172± 0.0060 


2746±140 


SE shell 


01:04:04.04 


-72:02:00.5 


Fig.E] 


12.8145± 0.0010 


0.1529± 0.0023 


3579±54 


2295l|3 



Velocity dispersion shift from the line centroid. 

Velocitytrtie is estimated velocity dispersion after accounting for the instrumental resolution. 



Table 4. Properties of Freshly Formed Dust in E0102 



Mofk^l 




''Dust Compositions (tem]")oratur(^ [K]. mass [M: ]) 


■''Mass (M. ) 


Total Mass (M ; ) 


A 


2.41(~398/165) 


MgSiOs (116, 3.80E-5), Si (150, 6.20E-5), C (55, 5.00E-3) 


^5.12x10-3 


" 0.015 


B 


1.98(~389/196) 


MgSiOa (117, 3.21E-5), Si (114, 4.4E-4), AI2O3 (=55, 2.06E-3) 


^2.53x10-^ 


^0.007 


C 


2.18 (~438/201) 


MgSiOa (108, 6.25E-5), Si (142, 9.21E-5), ^^203 (^U, 0.03) 


'^0.03 


<^0.09 


D 


2.16 (417/193) 


MgSiOs (124, 2.34E-5), Si (149, 6.3E-5), 










AI2O3 (60, 3.33eE-4), C (55, 4.350E-3), Mg2Si04 (100, l.lOE-5) 


^.80x10-3 


^0.014 



'^Ax^ is the reduced and "dof is the degrees of freedom. 

^Compositions in the best fit, where a few primary compositions are written in bold, and alternative dust compositions are in italics. A 
combination of compositions are in a regular font. 1 

^The errors are on order of 10-20%. 

'^The errors could not be constraint. 

^The errors are on order of 20-40%. 

^Mass is the mass within the IRS slit and total mass is the mass estimated over the entire SNR (see the text for details). 
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Table 5: Summary of Ne Column Density and Mass Estimates 



Ne ionization state 


from Excitation-rate equation 


from Shock model^ 




Model A" 


Model B'' 


Model 




[Ne I] (cm-2) 


2.4 X 10^^ 






1.5 X lO^s 


[Ne II] (cm-2) 


2.7 X 10^3 


4.4 X 10^"^ 


2.4 X 10^^ 


7.8 X lO^s 


[Ne III] (cm-2) 


6.8 X 10^2 


6.5 X 10^*^ 


1.5 X 10^^ 


1.2 X 10^^ 


[Ne IV] (cm-2) 


9.6 X 10^^ 






8.3 X 10^^ 


[Ne V] (cm-2) 


5.9 X 10^1 


6.0 X 10" 


6.0 X 10" 


1.1 X lO^^ 


[Ne I - V] total (cm-^) 


7.8 X 10^*^ 


1.1 X 10^5 


2.6 X 10^6 


4.0 X 10^'' 


Ne mass'^ (M©) in IRS slit 


0.22 


0.003 


0.07 


0.014 


total Ne mass^ (Mq) 


0.63 


0.01 


0.21 


0.042 



"Model A: ne=7 x lO"* cm-^ T = 612 K; Model B: ne=l x 10* cm-^ T = 350 K; Model C: ne = l x 10^ cm-^ T = 
215 K. See the text and Figure [T2l for details. 
'The best estimates (in bold) are from our shock model. 

'^Ne Mass is the mass within the IRS slit and total Ne mass is the mass estimated over the entire SNR. 



Table 6: Comparison of Elements detected'^ in E0102 and Cas A 



Wavelength 


SNR 





Ne 


Mg 


Si 


S 


Fe 


Ar 


Ca 


H 


He 


references" 


infrared 


E0102 


Y 


Y 




U 


U 


pN 


N 




N 


N 


this paper 




Cas A 


Y 


Y 








Y 


Y 








1, 2 


optical 


E0102 


Y 


Y 


Y 


pN 


pN 


pN 


N 


N 


N 


N 


3 




Cas A 


Y 


U 


U 


U 


Y 


N 


Y 


N 


Y* 




4, 5 


X-rays 


E0102 


Y 


Y 


Y 


Y 


U 


U 


N 


Y 






6 




Cas A 


Y 


Y 


Y 


Y 


Y 


Y 


Y 








7 



"Y: detected, U: uncertain, pN: probably not detected, N: not detected, * indicates that the emission is likely from 
circumstellar medium instead of from ejecta (see the text for details). 

''(1) Ennis et al. (2007), (2) Douvine et al. (2000), (3) Blair et al. (200), (4) Fesen et al. (2006), (5) Fesen et al. 
(1996), (6) Flanagan et al. (2004), (7) Hwang, Hoh & Petre (2000) 
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Table 7: Comparison between the Ejecta Masses of the SNR E0102 and Nucleosynthesis Yields 





Ne (M0) 


(Mq) 


Ar (M0) 


IR ejecta mass 


0.042 (from Shock model) 




0.78 






X-ray ejecta mass 


0.1 


- 2° (0.84*) 




0.3 - 6" 






total mass of IR and X-ray ejecta 


0.14 


- 2.4 (0.882'') 




1.08 - 6.78 






Nucleosynthesis Models'^ 


A 


B 


C 


A 


B 


C 


A 


13 M0 model 


0.028 


0.023 


0.064 


0.21 


0.151 


0.308 


0.0028 


15 M0 model 


0.039 


0.033 


0.049 


0.433 


0.368 


0.520 


0.0063 


20 M0 model 


0.257 


0.191 




1.480 


0.800 




0.004 


25 M0 model 


0.651 


0.594 


0.565 


3.000 


2.997 


1.386 


0.007 


35 M0 model 






1.366 






3.891 


0.040 


40 M0 model 




0.720 






9.110 







"From Flanagan et al. (2006). 

''Estimates from our shock model. Note that our shock model matches the infrared lines. 

' ^The model nucleos ynthesis yields are taken from [Thielemann et all l|l996f ) (A), iNomoto et al] (|l997l ') (B), and 
IChieffi &: Limongi I (|2004i ') (C). The best estimate and matching nucleosynthesis yields are written in bold. 
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Fig. 1. — Four ima ges of E0102 on the sa me ang ular scale (clockwise fr om upper left): in- 
frar ed MIPS 24 um (IStanimirovic et alJ bood) . HST (IFinkelstein et aD bood). ACTA radio at 6 
cm (|Amv fc Baillliad v and Chandra X-rays (IFlanagan et alJl2004l ). The brightest ring of radio 
emission traces the forward shocked material, and the brightest ring of X-ray emission traces the 
reverse-shocked ejecta. 



Fig. 2. — IRS slit coverage is superposed on a MIPS 24//m image (jStanimirovic et al.l 120061 ) of 
E0102. North is to the right and east is up. The large diffuse source located southwest of E0102 is 
the [H II] regions of N76. 
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E0102 




10 20 30 

Wavelength (/xm) 



40 



Fig. 3.— 



Spitzer IRS spectrum of E0102. 
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E0102 




10 20 30 

Wavelength (/xm) 



40 



Fig. 4. — The spectrum before (top), and the spectrum (bottom) after, local-background subtrac- 
tion. 



East 



NeV 24.3 micron 
OVI/Fell 26 micron 



West 



SHI 33.5 micron 

Siil 34.8 micron 
Ne iii 36.0 



^ '■ 



'I 



N76 ■ 

— sill 

Sis 



Fig. 5. — Merged LL position-wavelength images (left: LL2, right: LLl); X-axis is spatial direction 
along the slit covering E0102 and N76 and y-axis is wavelength. The image shows that E0102 emits 
continuum, [Ne V], [Ne II], and [O IV] lines. 
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Fig. 6. — Mosaicked three-color infrared images of E0102. Red: MIPS 24 (from Stanimirovic et 
al. 2006); Green: IRAC 8 /xm; Blue: IRAC 5.8 ^m images. The IRAC 8^m flux ranges from 0.15 
to 0.2 MJy sr~^. Possible detection of faint 8/um emission (seen in yellow) from the southeastern 
shell is suggested. The IRS slit coverage is marked as a green box. Note that the remnant is not 
detected at 5.8 nm. The b right emission at the bottom right of the image is from N76. X-ray 
contours (jGaetz et al.ll200d ) are superposed. 
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Fig. 7.— Line map images of E0102: (a) [Ne II] 12.8 ^um, (b) [Ne III] 15.55 /im, (c) [O IV] + [Fe II] 
26 /im, ( d) continuum 16.5 -20.5 //m, (e) continuum 27-31 /im, and (f) [Si II] 34.8 ;um. X-ray 
contours (IGaetz et al.l 120001 ) are superposed. The spatial resolution of [Ne II] image (from SL) is 
~2" and of the other images (from LL) is ~5.1". 
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Fig. 8. — (left) Spitzer [Ne II] intensity map (colors are from yellow, red and blue in strength of 
brightness) overlaid with HST [O III] contours, (right) Color images of X-ray (blue), optical [O III] 
(green) and infrared [Ne II] 12.8/im (red) emission. 




Fig. 9.— Mosaicked three-color images of (a) E0102 and N76, (b) E0102 only: [Ne II] (blue), 
16.5-20.5 nm continuum (green) and [O IV] (red). The image on the left is centered on R.A.= 
01^03™45.6^ and Dec.= -72°02'02.3", and the FOV is 5.2'. The image on the right is centered on 
E0102 and the FOV is 2.7'. 
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Fig. 10. — Line maps [Ne II] (12.8/im) showing blue- to red-shifted wavelength regions: 12.5-12.7 
/im (a), 12.7-12.8 /um (b), 12.8-12.9 /im (c), and 12.9-13.1 /xm (d). The maps are overlaid with 
the Chandra X-ray contours. The positions of the spectra in Figure [TT] are marked as boxes in the 
maps, (a), (c) and (d) for the locations of "E", "Scenter", and "center", respectively (see Table [3]). 
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Fig. 11. — Spectral line profiles of [Ne II] toward three positions (see Table [3]) in E0102: Center 
(blue), E (green), and Scenter (red). The line broadening in the Center spectrum is consistent with 
a velocity of 5226 km s~^ (see Table 3 for details). 
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Fig. 12. — Line diagnostic contours of the ratios of [Ne III] 15/35 /xm (red) and [Ne V] 14.3/24.3 
/im (blue). The observed ratios are marked with thick sohd hnes for the best value of the ratios 
[Ne III] = 54 and [Ne V] = 1.76. Note that the density is an electron density. The shaded region 
shows the range of temperatures and densities allowed by errors for the [Ne III] ratio, whereas 
the allowed physical conditions of [Ne V] lines are marked with the thick blue contour lines. The 
contours of [Ne III] 15/36 and [Ne V] 14.3/15.3 intersect at the solution with a temperature of 612 
K and a density of 4x10"^ cm~^ (marked as an asterisk). The column density of [Ne III] is marked 
for three different densities (see Table [5]) . 
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Fig. 13. — Energy diagram of [Ne V]. Five levels are shown to the right in a unit of eV, with a 
zoom-in of the three ground-stage levels to the left (with a unit of K) . 
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Fig. 14. — E0102 IRS dust spectrum superposed on a dust fit (Model A in Table H]). The continuum 
was fitted with dust compositions of Si, MgSiOs and Carbon. The compositions suggest that the 
dust forms around carbon-burning layers. The data and the total fit are shown in blue and thick 
red lines, respectively. Contributions from each type of dust are shown with the dust temperatures 
given in parentheses. The gas lines are also shown in the plot, but they are excluded in spectral 
fitting (see Figure [T5]l 
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Fig. 15. — The dust continuum superposed on the dust fit of Model B (Table S]) with the dust 
compositions of Si, MgSiOs, and AI2O3. The synchrotron continuum contribution estimated from 
the radio fluxes and spectral index is shown by the green dashed line. 
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Fig. 16. — The dust continuum superposed on a combination of dust models of Si, MgSiOs, AI2O3, 
Mg2Si04, and carbon. 
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Fig. 17. — Column density for each Ne ionization state. The column densities are derived from 
collisional ionization (squares) and from shock models (Filled dots). The low (5667 K, dash-dotted 
line) and high (9300 K, dashed line) temperature components are marked with the sum (solid line) 
of the two temperature components. The corresponding Ne mass is shown in the right side of 
y-axis. The horizontal lines (dash-three-dotted line) are [Ne I] and [Ne II] column densities and 
masses. 
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Fig. 18. — Ionizing radiation generated in the sh ock and incident on the gas in the photoionization 
zone for the case ~200 km s^^ and 3x the iBlair et al.l (2000) ram pressure. The spectrum 
is clearly line dominated and many of the lines are from O and Ne which have highly enhanced 
abundances in the ejecta. 
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Fig. 19. — Neon ionization (above) and gas tem perature (below) i n the shock for our model with 
shock speed 200 km s~^ and ram pressure 3x the lBlair et al.l (2000) assumed value. The ionization 
of Ne clearly shows the lag in ionization and recombination in the ionization (T > 10^ K) and 
cooling (10^ << r < 10^ K) zones. The region on the right, where T ~ 550 K, is photoionized 
and assumed to be in thermal and ionization equilibrium. The Ne^^ ion fraction is very low in the 
photoionized zone and so the observed [Ne V] 14.3 ^m and 24.3 ^m emission come from the cooling 
zone (an electron density of 10'^ — lO^cm"^) whereas most of the [Ne II] and [Ne III] emission comes 
from the photoionized zone. 
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Fig. 20. — Predicted [Ne V] line brightness and the ratio of 14.3/24.3 fim hnes depending on 
shock model parameters of shock velocity (indicated as numbers for the standard case runs), ram 
pressure (P ) and magnetific fields (B). Pfac and Bfac indicate the ratio of the value relative to the 



Blair et al.l (|200Cl ) values for the ram pressure and magnetic field. The observed Ne V lines are 



best fit by two shock models of {Pfac = 3, Bfac = 1) (blue line) and {Pfac = 1, Bfac = 0.5) 
(light blue line) with a shock velocity of 200 — 300 km s~^. The horizontal dotted line indicates 
the observed emission line ratio and the dashed lines show the formal 1 — a errors. The leftmost 
vertical lines indicate the 14.3 ^m surface brightness assuming that emission fills the field of view. 
If the unresolved emission regions only have a filling factor of 0.13 then the true surface brightness 
is that indicated by the vertical lines to the right. Our "best" model assumes that is the case. 
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Fig. 21. — (a) Nucleosynthesis yields of Ne mass depending on the progenitor mass, (b) Nucleosyn- 
thesis yields of Ar/Nc (dashed line), Ar/0 (dotted line) ratios and Ar mass (solid line) depending 
on the progenitor mass. 



